Even-odd behavior of conductance in monatomic sodium wires 
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With the aid of the Friedel sum rule, we perform first-principles calculations of conductances through 
monatomic Na wires, taking into account the sharp tip geometry and discrete atomic structure of 
electrodes. We find that conductances (G) depend on the number (L) of atoms in the wires; G is 
Go(= 2e 2 /h) for odd L, independent of the wire geometry, while G is generally smaller than Go 
and sensitive to the wire structure for even L. This even-odd behavior is attributed to the charge 
neutrality and resonant character due to the sharp tip structure. We suggest that similar even-odd 
behavior may appear in other monovalent atomic wires. 
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Atomic contacts have been generated by scanning tun- 
neling microscopes or mechanically controllable break- 
junctions 0]. It has been demonstrated that single 
atom contacts made of monovalent atoms such as Na, 
K, and Au have a strong tendency towards the quantized 
conductance Go (= 2e 2 /h). To explain the quantiza- 
tion, a ballistic transport model has been suggested |Q, 
where atomic contacts are modeled by jellium constric- 
tions, which are adiabatically connected to electrodes, 
analogous to quantum point contacts (QPC) in two- 
dimensional electron gas systems ||f|]. However, this 
analogy may be inappropriate for sharp tip structures 
of atomic contacts || , because the Fermi wavelengths in 
metallic systems are comparable to atomic spacings. Re- 
cently, Yeyati and his coworkers || have suggested a res- 
onant transport as an alternative explanation, based on 
their tight-binding calculations assuming the local charge 
neutrality. Despite several theoretical attempts |^-^] , the 
origin of the quantized conductance has not been clearly 
understood because the validity of simplifications such as 
the local charge neutrality assumption |(||| or a jellium 
electrode model || , which ignores the sharp tip geometry 
for electrodes, is rather unclear. 

In this Letter, to clarify the origin of the conduc- 
tance quantization, we calculate conductances through 
monatomic Na wires fl(| , taking into account the sharp 
tip shape and discrete atomic structure of electrodes. We 
perform real-space multigrid electronic structure calcula- 
tions within the local-density-functional approximation 
(LDA). Using the Friedel sum rule jll],[l2|, which relates 
the conductance to the density of states (DOS), we are 
able to reduce greatly computational demands. We find 
a resonant character in transport and a robust quanti- 
zation of conductances when the number (L) of atoms 
in wires is odd. For even L, conductances are not quan- 
tized, resulting in even-odd behavior in transport. This 
feature is very different from the length-independent con- 
ductances observed in ballistic quantum wires. 

Figure [l](a) shows a monatomic Na wire connected to 



two electrodes, which have an inversion symmetry. Each 
electrode is modeled by a cluster of M Na atoms in the 
bcc structure. The (111) direction of the bcc lattice is 
aligned with the wire with L Na atoms, which consists 
of 2 apex atoms in the clusters and L — 2 Na atoms be- 
tween the clusters. For the interatomic spacing d within 
the wire, we use the bond distance of do = 3.659 A in bcc 
bulk Na, and also test various different values for d. Using 
the inversion symmetry, the DOS can be decomposed into 
even (p e ) and odd parity (p a ) components. When only 
one eigenchannel contributes to transport |l],|9| , conduc- 
tances through the wire can be expressed in terms of p e 
and p using the Friedel sum rule |ll|,|l2[ ; 



G = — — sin 

h 



:(N e -N ) 
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where N e r \ [= J Ep dE p e ^(E)] denotes the number of 
electrons with even (odd) parity and Ep is the Fermi 
energy. With Eq. (|j), we can reduce greatly the compu- 
tational demand and analyze conductance characteristics 
in terms of the DOSs. 

We calculate self-consistently the electronic energies of 
the wire system using the real-space multigrid method 
JLlj within the LDA. Norm-conserving pseudopotentials 
are generated by the scheme of Troullier and Martins, 
and then transformed into the separable form of Klein- 
man and Bylander fTlj ] . We use a supercell geometry con- 
taining two electrodes, each of which has a cross section 
of 20.76 A wide. Electrodes in neighboring supercells are 
separated by more than 11.64 A, so that inter-supercell 
interactions are negligible. Using a grid spacing of 0.43 
A, we ensure that the total energy is converged to within 
10~ 7 Ry. To calculate the DOS, we use the Fermi-Dirac 
broadening of Et — 0.052 eV (600 K) for the occupation 
of each level. We find that the level splitting between 
even and odd parity states, which have negligible wire 
characteristics, is less than 0.003 eV. Since this value is 
much smaller than Et and the average level spacing near 
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Ep, we do not expect numerical errors in N e — N Q due 
to the charge sloshing problem [ jl6[ ] . 

The conductances calculated from the DOSs are plot- 
ted as a function of L for M = 64 and 95 in Fig. ||. In 
both cases G oscillates with L, exhibiting even-odd be- 
havior. For odd L, we obtain G = Gq since the difference 
(N e - N ) almost equals 1. For both M = 64 and 95, 
similar results are found, indicating that the finite size 
effect of clusters is negligible for odd L. For even L, on 
the other hand, N e — N a is not equal to 1 and conduc- 
tances deviate considerably from Go- Since the results 
for M = 64 and 95 are different from each other, the 
finite size effect seems to be more important for even L. 
At this point, it is beyond our computational capability 
to calculate conductances for larger cluster sizes. How- 
ever, based on the DOS analysis given below, we confirm 
that the qualitative feature of G being smaller than Go 
for even L is not an artifact due to the finite size effect. 

The DOS analysis provides useful information. From 
95-atom cluster calculations, the results for p e (E) and 
p (E) are drawn for the L = 4 and 5 wires in Fig. 0, 
and compared with p^(E) and p^(E) for their reference 
systems (hereafter denoted by the superscript R), which 
are defined by removing L — 2 central atoms in the wire 
without altering the inter-cluster distance. We find that 
p e (o) and pf( ) are almost identical over a wide range of 
energies except for near the Fermi level. From the pro- 
jected DOSs (PDOS) onto the L — 2 central atoms in 
the wire (here a sphere with the radius of 0.5 d is cho- 
sen for each atom), the difference between p e ^(E) and 
p\o) (-^0 ' s f° un d to be mainly caused by the existence of 
resonance states. Thus, unlike conventional QPCs, the 
monatomic wires have the resonant character, consistent 
with previous results . 

We estimate the spatial extension of the resonance 
states by comparing the charge densities of the wire sys- 
tem (L = 5, M = 64) with those for its reference system 
[see Fig. 0(b)] . Besides the major difference in the wire 
region, an oscillatory feature due to the charge rearrange- 
ment appears on a few layers in each electrode, in good 
agreement with the previous calculations (l7| . Thus, the 
resonance states are not strictly localized in the wire but 
extended somewhat into the electrodes. 

In the PDOS, there exists a half-filled resonance state 
at Ep for odd L (on-resonance) , while the Fermi level lies 
between two resonance states for even L (off-resonance) . 
The position of the resonance states affects N e — N a , and 
thus conductances. For demonstration, we use the rela- 
tion N e - N = ( N e - ATf) -(N - N*) , where A^ e R = N* 
because G R = in the reference system. Since each res- 
onance contains exactly one extra energy level p| , each 
filled resonance state with even (odd) parity increases 
N e ( ) - by 2 without altering N ( e ) ~ N f(e)> and 

thus increases (decreases) N e — N a by 2; the factor two 
denotes spin degeneracy. However, the filled resonance 



states are irrelevant because G does not change when 
N e —N a changes by ±2, as shown in Eq. ([!]). On the other 
hand, the half-filled resonance state changes A^ — N a by 
±1, and if there is such a state, G equals Go. When the 
Fermi level is located at the tails of resonances, \N e — N Q \ 
can be quite different from 1 and G < Go- This explains 
the correlation between the even-odd behavior of conduc- 
tances and the relative position of resonance states with 
respect to Ep. 

A recent tight-binding study || for L = 1 has demon- 
strated that the position of resonance states cannot be 
arbitrary when the local charge neutrality is assumed. 
Although this assumption is violated in atomic scale due 
to the charge oscillation in Fig. 0(b), we find that the 
charge neutrality does hold globally near the wire and 
constrains the resonance position for general L. Since 
Na is a monovalent atom, the wire system contains AA" 
(= L — 2) extra electrons, as compared to its reference 
system. On the other hand, we verify that the numbers of 
filled resonances in the two systems differ by AN/2 [fl9|| , 
which is a half- integer for odd L. This matching indicates 
that all extra electrons are in the resonance states near 
the wire. Thus, the alternation of on- and off-resonances 
with L results from the charge neutrality. 

The same even-odd behavior is expected for a wide 
class of wire structures, since the resonant character and 
the charge neutrality are common features. We exam- 
ine various wire structures such as stretched wires with 
d=l.l and 1.2 do and zigzag wires with the bond angles 
6 = 150° and 120° (d = d ) [see the inset in Fig. 1(a)]. In 
all cases, a half-filled resonance state is found for odd L 
and G = Go to within 2 %, while the off-resonance trans- 
port is realized for even L and G < Go- Thus, the same 
even-odd behavior holds for other wire structures. On 
a quantitative level, however, G is sensitive to the wire 
structure for even L. For the (L = 4, M = 64) system, for 
example, G/Gq is 0.61 for a linear unstretched wire, 0.50 
for a stretched wire with d — 1.2do, and 0.71 for a zigzag 
wire with = 120°. Recalling that in the off-resonance 
transport, G depends not only on the relative position of 
resonance states but also on the ratio T/AE, where AE 
and r denote the resonance spacing and width, respec- 
tively, nonuniversal G values for even L can be under- 
stood since the ratio depends on wire structure. For suf- 
ficiently large L, however, G becomes universal for even 
L as well since AE — ► and G — > Go as L — > oo [p0| . In 
this limit, the even-odd behavior disappears. 

Here we note that two key ingredients responsible for 
the even-odd behavior are the sharp tip structure and 
charge neutrality, which are believed to be common in 
other monovalent metallic wires. It is then very plausi- 
ble that the same even-odd behavior may appear in other 
monovalent metallic wires. Previous experimental stud- 
ies p],pi| showed that atomic contacts made of various 
monovalent metals have many properties in common. In 
agreement with our expectation, recent calculations |22] 
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for gold wires found that G is smaller than Go and sen- 
sitive to wire structure for L = 4. Experimentally, gold 
wires with L > 1 have been already reported |23|,^4| . Al- 
though these experiments did not exhibit an evidence for 
the even-odd behavior, it is premature to reject the pos- 
sibility of the even-odd behavior in gold wires because 
it is not clear whether uncontaminated gold wires with 
even and odd L are both produced in the experiments 
p5j. More experimental studies are required to clarify 
the existence of the even-odd behavior. 

We next discuss briefly atomic wires made of multi- 
valent atoms. Since multiple eigenchannels contribute 
to transport, Eq. (||) should be replaced by the general- 
ized formula Q, G = G J2 J J=1 sm 2 [(ir/2)(N e - JV )a 3 -], 

where J is the number of channels and X}/=i a j = !• 
One crucial difference from monovalent wires is that the 
new parameter aj does depend on wire structure. Thus, 
even when the charge neutrality fixes N e — N Q , G still de- 
pends on wire structure. This explains why multivalent 
metallic wires such as Al Q] do not exhibit clearly quan- 
tized conductances as in monovalent wires. We point 
out that for L = 1, the multichannel transport may be 
relevant even for sodium wires. We find that the cou- 
pling between two electrodes is not negligible even when 
the shared apex atom is removed. This implies that the 
direct transport between the electrodes, which is not me- 
diated by the common apex atom, is not negligible, in 
agreement with Ref. (9). For more precise conductance 
calculations for L = 1, the generalized formula should 
be used, which is beyond the scope of our paper. In re- 
cent experiments for Na atom contacts [^6|, the tail of 
the conductance histogram peak at Go was shown to be 
extended up to 1.2 Go, which may be interpreted as an 
indication of at least two additional eigenchannels |27]] 
with a small contribution of ~ 0.1 Go for each channel. 

In conductance calculations using finite-sized elec- 
trodes, the level broadening Et should satisfy the follow- 
ing constraints; to retain the metallic nature, Et should 
be larger than single particle level spacing Se (~ 0.02 
eV near Ep), while it should be smaller than T and AE 
(~ 0.3 eV). When Et varies between 0.03 and 0.15 eV, we 
verify that G for odd L is almost independent (to within 
a few percent). For even L, on the other hand, the reso- 
nance width r depends weakly on Et and thus the choice 
of Et affects G quantitatively although the even-odd be- 
havior remains robust. We also point out that T depends 
on the sharpness of the tip structure. For the pyramid- 
shape (100) tip, which is less sharper than the (111) tip, 
we find that T increases by about 40% and the overlap 
of resonances is enhanced, leading to weaker even-odd 
behavior and less accurate conductance quantization for 
odd L. Recent calculations || for an L = 3 wire with 
flat electrodes showed that the energy dependence of the 
transmission coefficient is almost negligible, implying an 
even larger T of about 1 eV. 



Finally, we note that when the exact inversion (or 
mirror reflection) symmetry is relaxed, the even-odd be- 
havior still occurs because the on(off)-resonance for odd 
(even) L arises from the sharp tip structure and charge 
neutrality. However, the conductance quantization for 
odd L is weakened, as illustrated clearly in Ref. ||. Re- 
cent experiments |2§| ] have demonstrated that the elec- 
trodes of atomic contacts tend to be aligned with sym- 
metric lattice axes. Here we remark that a similar con- 
ductance oscillation due to the resonance states is re- 
alized for heterogeneous systems j2(l, where wires and 
electrodes are made of different atomic elements. The 
origin of the resonance is however different; in heteroge- 
neous systems, the charge transfer between the wires and 
electrodes is likely to occur and the resulting Schottky- 
like potential barriers at both the wire ends generate the 
resonance states even if the electrodes are flat. 

In summary, a first-principles method is implemented 
to calculate conductances through Na wires, with the aid 
of the Friedel sum rule. Unlike conventional QPCs, the 
transport in monatomic Na wires has resonant character 
due to the sharp tip structure. Combined with the charge 
neutrality, the resonance states lead to robust quantized 
conductances when the number (L) of atoms in the wire is 
an odd number. For even L, the off-resonance transport 
is realized with G smaller than Go- We suggest that the 
same even-odd behavior of conductances may appear in 
other monovalent atomic wires. 
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FIG. 1. (a) The atomic structure for the L = 5 Na wire 
connected to electrodes, with the inversion symmetry. Each 
electrode is modeled by a cluster of M (=64) Na atoms, (b) 
Contour plot of the difference of the total charge densities 
between the (L — 5, M = 64) wire and its reference system: 
solid contours for 0.001 to 0.005 a.u. with the increment of 
0.001 a.u., whereas dotted and dashed contours for -0.0003 
to 0.0003 a.u. with the increment of 0.0002 a.u. Atomic 
positions in the wire are marked by gray dots. A zigzag wire 
is shown in the inset in (a). 
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FIG. 2. For unstretched linear wires with (a) M — 64 and 
(b) 95, G(L) (open circles) and N e — N a (filled squares) are 
plotted as a function of L. For G(l) marked by asterisks, its 
value may be modified due to the possibility of multi-channel 
transport (see text). 
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3. For a linear wire with d — do, pe(E) (solid lines) 



and po{E) (dotted lines) are compared with p^(E)[= p^(E)] 
(dashed lines). The Fermi level (Ef) is aligned at eV. Black 
circles denote the locations of the resonance states which are 
localized at the apex atoms of the reference system. Thick 
solid and dotted lines denote the projected DOSs (PDOSs) 
with even and odd parities, respectively. 
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